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Spin asymmetry of the ground states is studied for the trapped spin-degenerate (two-component) 
gases of the fermionic atoms with the repulsive interaction between different components, and, for 
large particle number, the asymmetric (collective ferromagnetic) states are shown to be stable be- 
cause it can be energetically favorable to increase the fermi energy of one component rather than the 
increase of the interaction energy between up-down components. We formulate the Thomas-Fermi 
equations and show the algebraic methods to solve them. From the Thomas-Fermi solutions, we 
find three kinds of ground states in finite system: 1) paramagnetic (spin-symmetric), 2) ferromag- 
netic (equilibrium) and 3) ferromagnetic (nonequilibrium) states. We show the density profiles and 
the critical atom numbers for these states obtained analytically, and, in ferromagnetic states, the 
spin-asymmetries are shown to occur in the central regions of the trapped gas, and grows up with 
increasing particle number. Based on the obtained results, we discuss the experimental conditions 
and current difficulties to realize the ferromagnetic states of the trapped atom gas, which should be 
' overcome. 
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Recent developments of laser trapping and cooling of atomic gases, which realized the Bose-Einstein condensates 
(BEC) for Alkali atoms Q, have opened up new interests for condensed atomic gases with multi-components: multi- 
• component BEC [^]|| and fermi degenerate (FD) gas and also bose-fermi mixed gas in trapped potentials 

One of interesting physics in such systems is the phase structures in the ground states. In general, multi-component 
systems show a variety of phase structures and phase transitions between them. A typical example is seen in the 
superfluid 3 He (a system with 2x2 components of orbital and spin angular momenta), where there exist more than 
'■^J ' two condensed phases resulted from different combinations of condensed components p"l| ] . 

In the condensed system of trapped atomic gases with multi-components, several studies have been done about 
new phases and phase transitions: the BCS states in attractive two-component fermi gas |l2|,[l3|], a transition to 
the superfluid states in the crossover region between BCS and BEC theories Ibjj-ITfl, phase separations of different 
components in the BEC |l^|l9[ , fermi gas j2Cj-|2^] and bose-fermi mixed system~ |23| , |24| . 

In this paper, we discuss the trapped fermi-degenerate atomic gas with two-components: m — ±1/2 (magnetic 
quantum number) states of spin-1/2 atoms, or two of the hyperfine states of fermionic atoms with larger total spin 
(F > 1/2). Thus, we denote these components as spin degrees of freedom (up and down components, or up and down 
atoms) in this paper. 

When the interaction between different spin components is weak, the stable ground state of the system has a 
symmetric state in spin-component distributions at T = where equal number of up and down atoms make similar 
fermi-degenerate distributions so that the total spin of the whole system vanishes (paramagnetic state). However, for 
sufficiently strong interactions between components, an asymmetric state with unequal numbers of up/down atoms 
can be stable because it can be energetically favorable to increase the fermi energy of one component rather than 
the increase of the interaction energy between up-down components. (The interaction energy between up- up/down- 
down components can be neglected due to the Pauli blocking effects.) We call it "collective ferromagnetic state" of 
two-component fermi gas. 

As the terms "paramagnetic" and "ferromagnetic" suggest, this theory is related with a mechanism of the metallic 
magnetism, originally proposed by Bloch and developed by Stoner et.al. poj ], where fermi particles are the electrons in 
conductive bands. It should be noted that, in the theory of metallic magnetism, a role of E- lnt is done by the Hartree- 
Fock exchange energy. Studies of "collective ferromagnetic states" of atom gas are also interesting as experimental 
testing grounds of theoretical ideas in magnetism. Recently, such ferromagnetism aroused new interests in neutron-star 
physics as a mechanism of the magnetar (neutron stars with strong magnetic fields) p6| . 

In the uniform system, such ferromagnetic states have been discussed for the atom gas of 6 Li in relation with the 
stability of the BCS states |0|. The ferromagnetic states have also been discussed on the trapped BEC where the 
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interactions between different spin components are the origin of the asymmetry p7] . In this paper, we treat the 
ferromagnetic states in the trapped and finite fermion system. 

In the next section of this paper, we formulate a set of equations for the two-component fer mi g as at T = and 
the stability condition of its solutions in the Thomas- Fermi (TF) approximations. In Section III, solutions of the 
TF-equations are analyzed algebraically, and a critical condition for the collective ferromagnetic states is obtained. In 
Section IV , we show the density distributions of the two-component fermi gas and discuss paramagnetic- ferromagnetic 
transitions in it. In Section [v], we summarize the results and discuss the experimental conditions to obtain the 
ferromagnetic states. 



II. THOMAS-FERMI EQUATIONS FOR TWO-COMPONENT FERMI GAS 

We consider a T = system of two-component Fermi gas trapped in an isotropic harmonic oscillator potential; 
densities of spin up/down components are denoted by p\ (r) and pi{r). To describe the ground-state behaviors of the 
system, we use the Thomas- Fermi approximation [ |12[]28|] , where the total energy of the system is a functional of the 
densities: 



d :i , 



r=1.2 



h 2 3 ( 

2m 5 



, n. 2 | 1 2 2 

wr ) 3 pix + -muj r p 



gpiP2 



(1) 



The m and uo in (yj) are the fermion mass and the oscillator frequency of the trapping potential. The last term 
in Eq. (|l|) corresponds to the interaction energy between different components of fermion, and the strength of the 
coupling constant g is given by g — AnT^a/m where a is the s-wave scattering length. The interactions between the 
same components are neglected; the elastic s-wave scattering is absent because of the Pauli blocking effects and the 
p-wave scattering is suppressed below ~ 100 pK. In the present paper, we discuss a system of repulsive interaction, 
so that the parameter g should be positive (g > 0). 

As we show later, the collective ferromagnetic states of two-component fermi gas show up in the case of large 
particle number (10 6 ~ 10 13 ), which validates the use of the Thomas- Fermi approximation. Also the validity of the 
approximation can be estimated from the smoothness of the mean-field potential: V e g — ^muj 2 r 2 + gp\.2- As a 
parameter that represent its smoothness, we can take a local de Broglie wave length A(r) = h/p(r), where p(r) is a 
TF local momentum defined by p(r) — y/2m(eF — V c g) (ejr : Fermi energy). As presented in |p9|] , a validity condition 
with A(r) is given by f(r) = |^| <C 1. Evaluating f(r) with the parameters given in the last section (for 40 K), we 
obtain f(r) < 10 -3 except the classical turning points. It also supports the validity of the TF approximations in the 
present case. 

To simplify Eq. (hi), we introduce the scaled dimensionless variables: 
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where £ = y/h/mui is the oscillator length. Using these variables, Eq. (|l|) becomes 
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The Thomas-Fermi equations for the densities ni,2 are derived from the variations of the total Energy E on rix^ 
with a constraint on the total particle number N: g^—(E — XN) = 0, where N is the scaled total particle number 
defined by 



N = Ni + N 2 = J2 J d3 



(4) 



The Lagrange multiplier A in the variational equation is for the fermion- number constraint, and related with the 
scaled chemical potential p through the relation pN — XN; using Eq. (0), we obtain 



p - 0.346 [-] A. 



(5) 
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It should be noted that, in Eq. (||), parameters (m, w, g) has been scaled out and no parameters are included except 
A. The Lagrange multiplier A is determined by the total fermion number N, so that N is the only parameter that 
determine the ground-state properties of the system. 

Using Eqs. (||) and (^) for the variational equation, we obtain the TF equations: 



n 'l + n 2 = A - 

The stability condition for solutions of Eq. 
functional: 



= M(x), nj + m = M(x). (6) 
can be derived from the second-order variations of the energy 
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Using Eq. 
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), we obtain the stability condition for the present case: 
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III. SOLUTIONS OF THOMAS-FERMI EQUATIONS 



In this section, we show solutions of coupled TF equations (^|) in algebraic form and discuss their stability based 
on the stability condition (||). 

Let's introduce variables s and t as ri\ = s 3 and n 2 = t 3 , then Eq. (0) becomes 



s 3 + t 2 = M, t 3 + s 2 = M. 
Alaking sum and difference of them, we obtain two equations equivalent with Eq. (^|) : 



s 3 + t 3 + s 2 + t 2 
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The factorized form of the second equation gives two alternatives: 1) s = t or 2) s 2 + st + t 2 — s — t = 0. 



In the case 1), s = t, the first equation in (|1C|) becomes s 3 + s 2 = M, which can be solved algebraically with the 
Caldano's formula. It includes only one positive root when M > 0: 
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For this solution, the stability condition M) gives s < 2/3, which leads to the constraint for M: 
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Thus, the stable symmetric solutions exist when < M < 20/27. 

Next, we take the case 2) : s 2 + st + 1 2 — s — t = 0. Now, this and the another equation in ( |Io|) are symmetric under 
the exchange between s and t, so that they can be represented by elementary symmetric polynomials, A = s + 1 and 
B = st: 



A 3 - A 2 - A + M = 0, A 2 - A - B = 0. 
The first equation has a real positive solution: 
, 9(M) 8 1 
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-108M + 44 + 12V81M 2 - 66M - 15 
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Substituting it in (|13| ), we obtain the solution for B: 
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The s and t can be recovered as two solutions of the equation: x 2 — Ax + B = 0. These solutions are generally 
asymmetric (s 7^ t). 

We can check the stability of the asymmetric solutions with Eq. (||), which gives a constraint for B: B = st — 
(nini) 1 / 3 < 4/9 . Then, from the second equation in (|l3|), we obtain that for A: 1 < A < 4/3. Differentiating the first 
equation in (113) by A, we find that M is a monotonically decreasing function of A within the interval 1 < A < 4/3. 
Combining these results, we obtain the stability range of M for the asymmetric solutions: 20/27 < M < 1. 

In summary, the stable solutions of Eq. ( [Tof ) are 

1. < M < § symmetric, 

2. |2 < M < 1 asymmetric, (16) 

3. 1 < M (no stable solutions). 

In case that 1 < M, Eq. ( [l0| ) have no stable solutions. In the TF variational equation -^-(E — AiV) = 0, we 
assumed the chemical equilibrium between fermion 1 and 2, and put their chemical potentials (Lagrange multipliers) 
equal: A = X 1 = X 2 . Accordingly, in 1 < M, we should say that the system has no stable solutions in "equilibrium". 
Thus, for 1 < M, we should take the nonequilibrium states where all fermions occupy one component (complete 
asymmetry). If we assume s — 0, t is obtained by solving the equation: 

t 2 = M, s = 0. (17) 



IV. DENSITY PROFILES OF TWO-COMPONENT FERMI GAS 



Using the results in the previous section, we can calculate the density profiles of the two-component fermi gas for 
any values of A. Corresponding to the classification by M — A — x 2 in (|l6|), we should divide A into three regions: 

20 20 

(i) < A < -, (ii) — < A < 1, (hi) 1 < A, (18) 

and discuss qualitative profiles of the density distributions 71^2 . 

(i) Paramagnetic ground states. In these cases, M(x) satisfies M(x) = X — x 2 < 20/27 for any value of x. Thus, from 
( |l6| ) , the density distributions are composed of the symmetric solution (|ll]) in all spatial region. Consequently, 
the ground states are paramagnetic in the region (i): n\{x) = TI2 [x). In Fig. |l](a), the density profiles are shown 
when A = 1/2. They monotonically decreases with the scaled radial distance x, and vanishes at the TF cut-off: 
xtf = vA = 1/a/2 in this case. These states occurs in the case of small fermion number. 

(ii) Ferromagnetic ground states in equilibrium. In these cases, we obtain M(x) < 20/27 in the outside region 
(x > y/X- 20/27), but 20/27 < M(x) < 1 in the inside region (x < yJX - 20/27). Correspondingly, the density 
distributions become symmetric in the outside region and asymmetric in the inside region, so that the ground 
states become ferromagnetic in the region (ii). Because the condition M(x) < 1 is satisfied, solutions are in 
equilibrium in all spatial regions, and densities are partially asymmetric in the inside region. In Fig. |l](b), the 
density profiles are shown when A = 4/5 for m (solid line) and ri2 (dotted line). The border between symmetric 
and asymmetric regions is given by xas = \/X — 2ti/¥l = y'8/135 in this case. In the outside (symmetric) 
region, the both lines overlap and vanishes at Xtf = 2/ v5- 

(iii) Ferromagnetic ground states in nonequilibrium. In these cases, density profiles are also asymmetric in the inside 
region (x < xas) as in the above region. However, in the most inside region (x < \A — 1), M(x) > 1 is 
satisfied and, according to the condition 3 in (|l6), the density profiles become nonequilibrium and are given by 
complete asymmetric solutions in (p"7|). In Fig. l|(c), we show the density profiles when A = 3/2 for m (solid 
line) and 722 (dotted line). In the most inside region (x < xeq = VX — 1 = l/y/2), the density profiles are 
completely asymmetric (n-i — 0), and, in xeq < x < xas = -\/41/54, they are partially asymmetric. In the 
outside (symmetric) region (a; > xas), two lines overlap and vanishes at x = xtf = \/3/2. Thus, the ground 
states are ferromagnetic and in nonequilibrium in the region (iii) which corresponds to the cases of large fermion 
number. 
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In Fig. |2j, we show the dependence of A on the fermion number N, which is calculated by Eq. (|4|) using the fermion 
densities 711,2 (x) obtained above; for the ground states (solid line) and the paramagnetic states (dotted line). From 
the above discussions, for A < 20/27, the ground states are paramagnetic, so that the both lines overlap. As can be 
read off in this figure, the critical value A = 20/27 corresponds to Nc = 0.53, which is the critical fermion number for 
the transition between paramagnetic and ferromagnetic states. 

In Fig. |](a), the variations of the scaled total energies E are shown against N, for the ground (solid line) and 
paramagnetic (dotted) states. The E can be obtained by Eq. ([?]) as a function of A using the fermion densities nx,i{x)\ 
combined with the iV-dependence of A (shown in Fig. |2|), we can obtain its iV-dependence. In N < Nc — 0.53, both 
lines overlap completely because the ground states become paramagnetic. In Fig. ^(b), we plot the energy difference 
between the ground and paramagnetic states. In this figure, we can find that, in N > Nc, the ground-state energy 
(solid line) shifts lower than that of the paramagnetic states (dotted line); it shows that the ferromagnetic ground 
states become more stable in this region. 

As can be seen in Fig. ||(a) when N > 0.53, the energy difference between ferromagnetic and paramagnetic states 
\E — -Eparal is very small in comparison with the total energy E: e.g. \E — E vajVa \/E ~ 0.01 at TV = 2.0. A large part 
of E consist of the stacked kinetic energy due to the fermi degeneracy of fermions. Roughly speaking, its size can 
be estimated from the total energy of the noninteracting fermi gas (g = 0) trapped in the same harmonic oscillator 
potential; i^non = ^^(^r^O 4 ^ 3 : which increases rapidly with increasing N. The E non is plotted also in Fig. |^(a) 
(dot-dashed line), the amounts of which are almost ~ 75% of E. 

However, to evaluate the scale of the energy difference, we should compare it with an one-particle excitation energy 
at the Fermi surface, which can be estimated from the scaled chemical potential jx defined in (||). In case of 40 K 
with the harmonic oscillator frequency u> = 1000 Hz (for other parameters, see the next section), Eq. (||) becomes 
fx r~j 10 _13 A. When N ~ 2 (A ~ 1 from Fig. ^), the ratio AE/jl becomes 10 13 . Thus, we should say that the energy 
difference between ferromagnetic and paramagnetic states are fairly large. 

In Fig.J|, we show the fermion number asymmetry (Ni — N^/N = (iVi — A r 2)/(Ai + iVjj) against the total fermion 
number N. 



V. SUMMARY AND DISCUSSIONS 



We discussed the possibility of transition to ferromagnetic states in the two-component fermi gas using the Thomas- 
Fermi approximation from theoretical points of view. Based on the results that we obtained, let's discuss experimental 
conditions and also difficulties to observe ferromagnetic ground states in the trapped atomic gas. We hope that these 
difficulties are overcome in future developments in experimental technics. 

As shown in the previous section, the ferromagnetic ground states become stable when N > 0.53; using Eq. (Q), 
the unsealed critical number Nc becomes 

0.53 /A 6 1 _ (i " 



^~ 0-346 =L5 UJ ' (19) 

As an example, we take the 40 K atoms (mass m = 0.649 x 10 _25 kg) trapped in the harmonic oscillator potential 
with u = 1000 Hz. For the scattering length, we take the value a = 169a^ (fiB- Bohr radius) given in |30|]. Using 
these parameters, we obtain Nc ~ 10 13 . In recent experiments, the trapped Fermi-degenerate gas has been performed 
up to ~ 10 6 atoms, so that the larger trapping potential are necessary for realization of the ferromagnetic ground 
state than currently used one. In addition, because of the high central density (~ 10 17 cm~ 3 ) for ~ 10 13 atoms, the 
inelastic/multi-body scattering processes in them becomes important, which might destroy the trapped atoms before 
they reach the required density. 

There can be several possibilities for the reduction of Nc- For example, if the scattering length can be increased 
by the Feshbach resonance for 40 K, which has been observed experimentally the value of Nc decreases and 
the ferromagnetic ground states can be obtained in small fermion number: simple estimation gives, a = 820a^ for 
Nc ~ 10 9 and a = 2600as for Nc ~ 10 6 . However, in current experiments, the tuning into Feshbach resonances is 
done by applying the magnetic field, which should make the energy difference between the up and down states and 
makes the ferromagnetic transition into a crossover in the spin asymmetry. The ways to tune Feshbach resonances by 
non-magnetic external forces (e.g. electric fields) are preferred for the present purpose. 

The use of heavier elements, e.g. Sr or Yb, is also effective for the ferromagnetic states. We hope that the 
combination of these methods may lead to the experimental achievements. 

We also comment on the process to perform the ferromagnetic states. There exist two possibilities: 
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a) The experiment starts with the magnetic held in some direction which makes the spin asymmetry (e.g. Ni > N2). 
Then, the magnetic field is switched off adiabatically in the cooling process. When the number of the remaining 
atoms is enough large, it shows the ferromagnetic states. 

b) The experiment starts with the symmetric trap (N\ — N2), and, after some relaxation time elapses, the atoms 
release their spin angular momenta and become the ferromagnetic state. 

The case a) is similar with a standard processe for observing the phase transition in the ferromagnetic materials. In 
case b), the spin relaxation time is considered to be the same order with that of clusterization, so that it might be 
difficult to observe the ferromagnetic transition within the lifetime of the atomic gas. 

Finally, we comment on the spatially phase-separated states in the trapped fermionic gas in the case of the large 
particle number and interaction strength ptj|-2^ ]. They correspond to the phase separation in the uniform system 



(two-phase coexistent region), discussed in [|12 . The naive TF calculations give smaller energy for the ferromagnetic 
states, but the energy difference is very small. The competition/coexistence of these states should be an interesting 
problem in future. 

We are very grateful to T. Maruyama, T. Suzuki and T. Miyakwa for fruitful discussions. Special thanks are also 
for T. Tatsumi for introducing us the physics of collective ferromagnetism. 
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FIG. 1. Density profiles of fermions, n\ (solid line) and 712 (dotted line) for (a) A = 0.5, (b) 0.8 and (c) 1.5. The scaled 
densities 711,2 an d the scaled radial distance x are dimensionless and defined by Eq. (^). The A is a (dimensionless) Lagrange 
multiplier. The (a)-(c) correspond to a paramagnetic state (i), a ferromagnetic state in equilibrium (ii), a ferromagnetic state 
in nonequilibrium (iii) each other. The scaled fermion number N\,2 are N\ — N2 = 0.089 for (a), (Ni,Ns) = (0.33,0.31) for 
(b), and (3.50,0.68) for (c). 

FIG. 2. Lagrange multiplier A against the scaled fermion number N. The solid and dotted lines are for the ground and 
paramagnetic states. The A is introduced as a Lagrange multiplier for fermion-number constraint, and the N is denned by 
Eq. (Q). Both quantities are dimensionless. 

FIG. 3. (a) The scaled total energy E against the scaled total fermion number N — Ni + N%: the solid line is for the 
ground states and the dotted one for the paramagnetic states. The dash-dotted line is for the total energy of the noninteracting 
fermionic system, E non = ^7r 2 (^-iV) 4 ' 3 . (b) The energy differences from that of the paramagnetic states, E — E paT& . The 
solid and dotted lines are for the ground and paramagnetic states. 

FIG. 4. Variation of the fermion number asymmetry (jVi — N2)/N against the scaled total fermion number N. The solid 
and dotted lines are for the ground and paramagnetic states each other. 
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